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1. INTRODUCTION 


One of the important problems of molecular chemistry concerns the evalua- 
tion of the correlation energy of molecules and the study, in general, of cor- 
telation effects in many-electron systems. The motions o/ the electrons in 
atomic or molecular systems are not independent of each other but are 
dependent on (or correlated with) the positions and spins of the other elec- 
trons. This correlation now is of two kinds—the first one arising from the 
limitations imposed by the Pauli Exclusion Principle which forbids two 
electrons of the same spin to stay in the same state, and the second kind 
arising from the strong Coulomb repulsion experienced by any two electrons 
when they try to approach each other closely. 


The standard method of solving many-electron problems in quantum 
mechanics is by means of the Hartree-Fock equations, which are the varia- 
tional equations of the Hamiltonian operator for a wave function that is a 
determinantal expression in the one-electron orbitals. Since the wave 
function is a determinant in the one-electron orbitals of the different electrons 
of the system, it will vanish when the co-ordinates of any two electrons 
having the same spin become identical and thus the Pauli Principle is implicit 
in this scheme. Or in other words, the Hartree-Fock equations suggest 
and bring within their scheme the correlation between electrons of the same 
spin, but because of the o~e-electron approximation they fail adequately to 
take into account of the Coulomb repulsion effects between the electrons. 
Experience has in fact shown that the H. F. equations give a much higher 
value for the energy of a molecular system than the experimentally observed 
one; this is to be expected because the H.F. equations neglect the Coulomb 
correlation between the electrons and would consequently contribute a much 
higher value for the repulsion energy of the system than is actually possible. 


Coulomb correlation between the electrons can be brought into the 
theoretical formalism in several ways, but the simplest and most direct method 
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of them all is to multiply the wave function by a correlation factor g (r,2, r,s) 
...+) which is symmetric in the inter-electronic distances ry., r3...... » ete. 
Such a method has been found to be very effective for two electron systems 
like He*ions and the hydrogen molecule, and it has been shown by Hylleraas, 
Léwdin and Redei that by using simple correlation factors of the type es» 
or (1 + ar,.), the energy for the Helium atom could be improved much 
beyond the value yielded by the H.F. equations. The introduction of the 
correlation factor modifies the field in which an electron is moving, and 
each electron moves in the average potential field of the remaining electrons 
subject to the condition that no two electrons can approach each other 
closely. The wave functions of the different electrons moving in such-a 
correlated field will be different from the one electron orbitals obtained by 
solving the H.F. equations, and an important problem in the study of cor- 
relation effects is to determine the one-electron orbitals in this case. “The 
equations determining these orbitals are given in Section 4 by making use 
of a formula given by Léwdin; in Section 8 we give the correlated Hartree- 
Fock equations for non-stationary systems. In Section 4, we have’ given 
the integro-differential equations satisfied by the generalised density matrices, 
and from these it is shown that one can obtain an expression for the-energy 
matrix of the system which will be useful in determining a correlated Thomas- 
Fermi destribution. In Section 7, some remarks are made on the effects 
of correlation by a study of the general equations for the simple case of two 
electron systems. 


2. DENSITY MATRICES 


Let ¥ (x,, Xo, ...,Xy) be the normalised eigenfunction of the system 
so that we have 


f P*¥dr = 1 1) 
where 


f dr = J dx,dx, eee dxx 


denotes integration over all space co-ordinates and summation over: spin 
co-ordinates. The generalised density matrices have been defined. by 
Léwdin as follows: 


Pr (x,' | x,) 


x Nf P* (x,'xXq ... Xn) F (xX, ... ¥n) dey ..: Dey 
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P’(x,'Xq | X12) 


= (3) f P* (x,'Xq'X3 ... Xn) PF (erty... Hy) dy... dey 


T'(x,'Xq) ... Xp’ | 4X2... Xp) ; 
-(*) f W* (x1'Xq 2. Xp'Xpia es Xn) WY (%y.-. Xw) 


X AXpis --- Un. (2) 
For a system of N electrons, ¥ is an anti-symmetric function of the co- 
ordinates (space as well as spin) of the electrons and we choose.for Y the 
following determinant: 
$y, (x%y) de (x) .-. by (%1) 
$y (X2) po (X2) ... ty (Xe) 


GO OOOO: 0. © 0 Oy OOo 01:6: O80: OO @ C181 000 


$y (Xn) te (XN) -.- Py (Xn) 


Y = (N!)?* (3) 


We assume for the sake of simplicity that the ¥’s form an ortho-normal 
system so that we have 


f da* (x) de (x) dx = dix. ie: (4) 
Now 
wey — L I det || p;; 
= NTP > NY EH PG (5) 
where 


N 
pi = z= by * (xi) dK (Xj). 
For future applications, we need also the quantity p’;; defined by 


pig = z be * (xi') ke (xj). : (6) 


ke1 
Since the ¥’s form an ortho-normal set, we have 


J p'inpnjaxn = p’ij- , - ; (7) 
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The integral of p over the co-ordinates of the various electrons has been 
given in Mott and Sneddon.> By adopting a similar procedure for the 
evaluation of the integrals and making use of (7), we can show that | 


P(x,’ | x) = p’'un = py 











ti. gf 
P (x4'x2! | Xy%2) = ui 9 a _ 1 Pa 
Pu +++ Pup 
P’(x'%q' «+ Xp’ |X... Xp) = ry vette eens = air (8) 
P'p1 - ++ P’pp 
where we have written p’p for the determinant in the above line containing 


p tows and p co'umns. 


There are several problems that deal with time-dependent processes,’ 
and the density matrices in such cases can be defined analogous to (2). Thus 
if Y (x,x_ ... Xyt) is the normalised time-dependent wave function for the 
system, we define the density matrix of order p as follows:— 


TD (%1'Xq! «.. Xp't! | XyXq ... Xpt) 


ion e P* (x4'Xq) ... Xp Xp... Xnl’) FP (xy... xt) 
Dp AXp+1 eee dxx. (2 a) 


If as before we choose for Y a determinant built from one-electron orbitals 
so that 


Y (x, ... Xyt) = (N!)-? det {y, (x12) ... by (xyt)} (3 a) 
where the orbitals y,, ¥., ...,% , satisfy the orthonormality conditions 

J de® (xt) Pe* (xt) dx = di (4a) 
then it follows that 





, 


I (x,'x3 eee Xp't’ | xy eee Xpt) = es . (8 a) 


In the above 


Pu P'xg --- P 1p 
Pp’ P'az.--> 





ng 


Ss" 


1e 


| 
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and 
pig = p (xi't'; xjt) = ZT pe* (x4'0') de (jt) (6a 
: 


Further if G(x,x_ ... x1) is a symmetric function of the co-ordinates o 
l electrons, it can be shown (see Appendix I for proof) that 


p'n 
NI f G (x,%2 --- x1) ys dx, ... dx 


l 2P'Lin (N — /— 1) 
= nl f G (x, eee x1) “— dx, oe dX 1.4 on a - 
P11 a ee Pil p 1141 
GG... ssi Bang Gn aie caael on 
x dx, ... AXta, Pile osteo se ca Pll Pile |: (9) 
2p L+3)1 ptt g 
rT; pot —" 








3. THE ENERGY OF THE SYSTEM 


The Hartree-Fock equations are the Euler differential equations of the 
variational principle 


8 f P* H¥dr = 0. (10) 


In view of the one-electron approximation, the Hartree-Fock equations 
do not adequately take into account of the Coulomb correlation, and we shall, 
following Léwdin, take the correlated wave function of the system as the 
product of a correlation factor which is a symmetric function of the inter- 
electronic distances and the Slater determinant. The wave function ® for 
the system then becomes 


® = g (rig, Ts...) ¥. (11) 
For g, we shall choose the form 
S(nans--J=atB sz fry) (11 a) 
i<j 


where the functions /(rjj) are supposed to be known. The simple functions 
ry, or (e47s: — 1) may be regarded as good choices for the correlation function 
f(r). One of the constants, say a, can be determined from the condition 


f ®* Odr = | (11 5) 
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and B can then be chosen as the best parameter minimising the energy. Now 
the best one-electron wave functions yj; (x) (i = 1, 2, ..., N) can be determined 
from the variational principle 


8W = 5 f O* HOdr = 8f ¥* H.Vdr = 0 


(12) 


(13) 


where 
H,. = gHg. 
Before writing down the variational equations we shall first fix the form 
of He. 
Now 


He = ‘ +P Zz reo} 1 Ej + 4 2 F (xj, “| 


i<j 


x jee = re 


i<j 


— 12 Ei +4 5’ FQ no| 


+ of = [f(rst), EW] 


s=1 «8,8 
s<t 


+208 5 ZF frst) F (xi, Xx) 


st ik 
e<t t<k 


+e SES feE WS(se) 


st ik a’t’ 
s<t i<k a’ <i’ 


+ 2 EF ZF frst) F Xi xn) f | 


= 7,+ T.+ T; + T,. 


(14) 


In the above [u, v] denotes the operator (uv + vu). Consider now the term T;. 
The summation here contains N*(N — 1)/2 terms. Of these, the number 
of terms for which i, s and ¢ are different is 3 (¥) while the number of terms 


‘nvolving two indices only is N(N — 1), 








Now 
lined 


(12) 


(13) 


Tm 
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Let us write 
T, (i, 5) = oB [f(ris), EC) + E(S)] 
Te (i, s, t) = oB {[f (rst), E@) + [fCrit), EO] 


+ [f (ris), E (O}- (15 a) 
Then 
T2 = a T (i, s) + 2 T, (i, S, t). (15 5) 


Consider next T;. The summation over i, k, (i< k) and‘s, t(s < ft) 
contains N?(N — 1)?/4 terms. Of these, there are (4) (%) terms of the 
type rst F (xj, xx) involving four different indices, 6(%) terms of the type 
F (xj, Xt) ris involving three different indices and (3) terms of the type 
rik F (xi, Xk) involving two indices only. If we write 


Ts (i, 5) = 28 f (ris) F (xi, Xs) 
Ts (i, 5, t) — 2a8 [F (xi, xt) (f (ris) +f (rst)) + F (xi, Xs) 


x (frit) + F (rst) + F Os, x0) (Ss) + Fd) 
T; (i, k, s, 1) = 208 [ f (riz) F (xs, xt) + five similar terms arising 
from the permutation of i, k, s and ¢], 
then 
T; = p> T:459)+ SF T3(is,t)+ 2 Ts3(i,k,s,t). (16) 


ik,s 4,k,8,t 
t<k i<k<s i<k<s<t 


We shall finally consider T,. The first term of T,, involving the operators 
E;, consists of N?(N — 1)2/4 terms and these can be split up into 30 (%) 
terms containing five different indices, 48 (Y) terms containing four differ- 
ent indices, 21 (*) terms containing three different indices and 2 (3) terms 
containing two different indices. Similarly the second term of T, involving 
the operator; F (xj, xx) consists of N?(N — 1)%/8 terms <nd the terms invol- 
ving 2, 3, ..., 6 different indices in this are given by the relation 


N= 99%) + 10M) + 14%) +24(3) 


+ iy] (17) 
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Let us denote by T, (i, 5), T, (i, 5, t), Ty (i,k, 5, t) ... Ty (i,k, 5, t, s’, t’) the 
terms in T, involving 2,3, ...,6 different indices. Then 













r= FUG)+..+ FEE UGK Sts,0). (18) 


9 st =a’? 
i<s iccectevcy 


We shall now write 
H,! = a* E(i) 
He? (i, s) = a? F (xj, xs) + Te (i, s) + Ts (i, 8) + Ty (i, 5) 
H,? (i, s, 1) = T, (i, s, t) + Ts (i, 5, t) + Ta (i, 5, 0) 













H,° (i, k, s, t, s, t') = T, (i, k, s, t, s’, t’). 






Then obviously 


He= J) He (+5 1 Dy Be (i,s) +. 
‘- 2) 2) BEG. . (20) 


We may note that the functions He! (x,, x3, ..-, x1) are symmetric in their 
arguments. 












Now 
W = f &* Hd; = f P* Headz. 
From (20) and (8) we have 





(21) 








W = J He! (xy) P(x’ | x) dx, + § He? (x1%2) PF (x4'x's | x1%9) dx dx, 
oe Hg? Ggtg. . ste) PGay .. 0g" | Xe... te ay... ite 
(22) 


where we are following the convention that, in the integrands, the operators 
He! (x1), He? (%1, X2) ... operate only on the unprimed co-ordinates x,, ..., X¢ 
and that, after these operations have been carried out, we have to put 
Xy' = X13 X_' = Xe; ...3 X—’ = Xg before the integrations. 







We see from (22) that the energy of the system is a functional of the 
density matrices alone, 
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he 4. THE CORRELATED HARTREE-FOCK EQUATIONS 


Léwdin has determined the one-electron functions 4,, #., ...,%x 
8) that make the determinant Y the ‘ best’ approximation for the eigenvalue 
problem 


Qop. ¥ = WP (23) 
where Qop, is a many-particle operator expressible as 
Oe 1 i 1 é 
Qop. = Ro + 2Q3 + 71 245 + 31 > Qijk +... . (24) 
‘ éj 4jk 


The prime in the above indicates that the summation excludes terms having 
9) two or more indices equal. 


The equations determining the functions ¥,, %.,...,%, are given by 
n(x) He (Xa) 
Qype (x) + f Qi. rt P Xe 
p(X’, X1) pp (X2'Xa) 








tie (x) =e (Xe) ~~ He (Xs) 
0) + 71 f Dyas | p(%a'X1) pp (X2'X_) pp (Xa'X3) | Axadxs 
p(X3'X1) pp (Xs'X2) ~—p- (3's) 





N 
=2 Ait (1). (25) 
Now H, is a many-particle operator similar to (24) and thus the best one- 
) electron orbitals that minimise the energy of the system can be obtained 
from (25) by replacing Q2op. by He. By expanding the determinants in (25) 
along the first row and denoting by p’);,; the co-factor of p’,; in p’;, we can 
write the integro-differential equations determining the correlated orbitals as 
pa 1 
) He? (x1) xk (x1) + d—1)! f He! (x,%2 ... *1) 
I=2 
" : 
. x ) P't; ke (xj) dxy ... dx 


N 
; — DF reiti @e) = 0 (k =1,2,...,N). (26) 


(a1 
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The conjugate complex of the above equation is given by 


He! (x1) oe * (x4) + a (1 . i)! f H,! Se 03 
l=2 


x > 01: 5b * ke (x4) Axe... dx 


j=1 


N ‘ 
“ > Anda* (x) = 0 (27) 


Equations (26) determine the correlated orbitals 4; (x;). They are 
a generalisation of the Hartree-Fock equations and since they explicitly take 
into account of correlation effects, they can be expected to give better theo- 
retical values for physical parameters like energy, than the simple Hartree- 
Fock system of equations. 


From (26) we have 


_ P a a7 f $i * (xy) He! (yy .., 2 


x s P'1; je (x5) dx, ... dxy (28) 


and from (27) we see that 
Aki = A* ix. (29) 


As in the H.F. scheme we see that the Ay; are a set of arbitrary constants 
forming a Hermitean matrix. Again as in the H.F. scheme, the set of orbitals 
ty, to, ..., Hy are not uniquely determined as the equations (26) are in- 
variant under an orthogonal transformation of these functions. The ele- 
ments p’j; of the density matrices (and consequently the density matrices 
themselves) are, however, uniquely determined. f 


Multiplying (26) by %,*(x,’), summing over k and remembering that 


Ul 


1 z — ’ ’ 
7 P'1:1jP' 17 = (xy, gee SI leave (30) 


j=1 














7) 


ire 
ke 


O- 
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we get from (26), (28) and (30) that 


6 


S Lf He (x, .<- TG Ww, .. alm. eee... a 


te1 


= Pa (I aa f PurHe! (XyX2 eos x1) 


Pip Piz +++ Pal 


4) OS ORs Ax dx... dx\. (31) 


Plu Ple-+- PU 


Multiplying (26) by 4,*(x,) and subtracting from it the corresponding 
equation for the conjugate of 4, (x_) multiplied by 4%, (x,) and summing 
over k, one gets 


De cee 


j=1 


Pai Paz cee Pal 


Por Pa2z---+ Pal dx, ... ax; 


eee eee eee eee ears 


_ Bi a i)! f H.! (x, ... x1) 


t=1 


Pia P1n--- Pil 


x Pea P22--- Pal dx, ... dx (32) 


Pla Ple--- PU 


Equations (31) and (32) are integro-differential equations governing the density 
matrices. When correlation effects are neglected, we have g = 1 and there- 
fore a = 1 and 8 = 0. It can easily be verified that in this case equations 
(26) reduce to the ordinary Hartree-Fock equations. Further, if we write 
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F (x, X2) = e* | ryg and B(x,) = J F (x, x2) p (xa, X2) dxq, then equation (32) 
reduces to the well-known equation.* 

{E (x) + B(x) — E(x’) — B (x’)} p (x, x’) 


— J [F (x, x”) — F (2’, x”)] p (x, x”) p (x’, x’) dx” = 0. (33) 


If we write 


a 2 
E(x) = om + V(x) = — ge Va? + V3), 


equation (31) becomes 
h? N\e 2 , , 
~ gain J 8 — ¥) Vato (a, x) de’ + (V + B)p (ex) 


= J F (x, x’) | p(%, x’) |? dx’ + Jp’, x) E(x’) p(x, x’) de’ 
+ J F (%’, x") | p(x, x’) |? p (%", x”) dx’dx” 
— SFO’, x") p @, x) p (x, x”) p (&’, x’) dx'dx’. (34) 
The above equation can be used to determine the charge density p (x, x) 
of the electron cloud of the atom. 
5. THE ENERGY MATRIX 
Let us now consider in greater detail the left-hand side of equation (32). 


We have 

6 , 
LHS. = He @eat D) gay J 
t=2 


Pai Pan --+ Pal 


rabicah dx, ... dx}. 


x eae + 








3) 
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Now since H,! (x, ... x1) is a symmetric function of x,, X2, .-.,X1, it can 
be seen by simple transpositions (x2, xx) (k = 3, ...,/) of the integration 
variables x2, Xs, ..., Xj that the (/ — 1) terms in the second summation of 


the above integral are all equal to f He! (x, ... %1) p’aah’l; 122X2 . -- AX1- 
Thus, 
P l 
LHS. = He! () 'a + gi f He! (x, ... xD 


x {p’a1P’1; u +@U— 1) P' a2P'l; 13} dx, ... ax}. (36) 


Now pij = p (xi, Xj) = Z du* (xi) oe (xj) can be regarded as the (xi, xj)-th 
element of a continuous matrix p. We shall define a matrix K (x9, x) 
by means of the following relation: 


K (2, X,) = K* (X2, 1) + EA re i! K! (x, x;) (37) 
Is2 

where 

K* (x2, %1) = J 8 (x, — 9) He! (, Pr) 8 (4 — x) dy 

K! (x9, X;) = 5 (X2 — x1) B! (x) +@¢—)) J He! (x ... x) 

x Pl: 120X3 — (38) 

and 

B! (x,) = f He! (x1... 1) p'ty1n Xn. Aen. (39) 
We shall define the (al)-th element of the product by pK! by 

(pK!)a, = JK! (x:, X1) P (XasX2) dX, (40) 


The definition is the same as the usual matrix multiplication law, but we take 
care to write p(x,, Xx.) after the operator K (x2, x,) so that the latter can 
operate on p(x,, Xz). We can then show the left-hand side of equation 
(32) to be the (al)-th element of pK. We have in fact 


(pK) a1 — j 5 (x — Xs) BI (x;) Pp (Xe Xe) dx» 
+ (1 — 1) J He! (xy... 1) p'tz 12'ag -- AXg ... dy 


t 
= B! (x,) par + J He! (x1... x) ES pts np akdX, ... dx 
k=2 
Pat Pas eee Pal 


ei TO iii | 
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Further He' (x;) p’a; is obviously the (a!)-th element of pK!. Thus the left- 
hand side of equation (32) is the (a’)-th element of pK. Similarly 


l 
R.H.S. = He! (x,) p’ + om + aa 


Pia Piz... py |* 


x Pe P2--- 


P'lg--.P 


= He! (xq) Par’ Sa d-— Dp! 


X JP (xpxe ... Xn) He! (xgXe ... x1) 
x Fez... eae, .... dite 


= He’ (Xq) Par’ + Sao Dt 4 1)! 


X J P* (xexs «.. Xy) He! (x—%q ... X 
xX P (xyre ... Xy) dxg ... dxy 
since H,! is a real operator. 
Thus we have 


, l 
mas meant Doty era 


P'ar Par -++ Pal 


x P'an Pag +++ Pal 


, 


Plao--s PU 


(Kp)q: = 2 (1 r. 1)! f K! Gx, X2) p’ (Xa, x) dx. 


> 1 ee 
= SF ty BO n+ DSH 


* P'L: gaP 2:4X2 eee dx}} 
= R.HS. 
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as can be seen-by expanding the determinants in (43) along their first columns. 
Thus the equation (32) is the (al)-th element of Kp — pK. Hence we can 
write 

Kp — pK = 0. (45) 


The second term in Ky.(x2,x,) subtracts from the matrix B-(x2, x,) = 
§ (x, — X,) B (x,) the physically irrelevant terms corresponding to the action of 
an electron upon itself and at the same time accounts for the exchange effect. 


We can in fact express the energy of the system in terms of the matrices 
K! (x, X2) we have 


l 
(29) - 


where D denotes the diagonal sum or spur. 


The matrix K is a generalisation of the matrix (E + B — A) (see Frenkel, 
pages 428-36) for the case when correlation is introduced into the theoretical 
formalism, and reduces to the latter when correlation is neglected. It can 
thus be regarded as the energy matrix for the system in analogy with the 
theory of the density matrix based on the Hartree-Fock equations. The 
importance of the energy matrix lies in the facility with which it enables one 
to pass on to a representation in the phase space and thus to obtain a semi- 
classical expression for the density of the charge cloud of the system. It is 
well known that by. transforming the matrix K (x, x2) into one {K (x, p))} 
involving the position and momentum of the particles, one can obtain 
the so-called Thomas-Fermi-Dirac equation that includes exchange effects 
besides. Thus the matrix K (x, x.) defined in (37) can be used asa con- 
venient starting-point to derive the equation governing the charge cloud of a 
molecular system; the details of this transformation theory leading to a 
correlated Thomas-Fermi charge distribution will be published’ separately 
in a different paper. 


6. A GENERAL CORRELATION FACTOR 


In the above discussion, we have chosen a correlation factor which is 
a linear function of / (r;;) where rjj is the inter-electronic distance between 
the i-th and j-th electrons; each term of g thus contains the co-ordinates of 
two electrons only. While this type of correlation factor is the simplest 
to deal with and can certainly be expected to improve the H.F. equations, 
it does not take into account of the multiple correlations connecting the 
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positions of different electrons at the same time. A general correlation 
function is either a polynomial in the (¥) inter-electronic distances or a 
convergent power series in them. Let us suppose that H, (= gHg) is of the 


form 
He = » px He! (uy a oni 47 


The correlated H.F. equations in this case become 


N t 
l 
> ay f He! (Xu) Xp «+ +5 20) > P'1; ike (5) dxe ... dx 
N 
- S Anis (x) = 0 


The difference between (48) and (26) consists only in the summation 
for /; while in (26) the summation for / runs from | to 6, in (48), it runs from 
1 to N. The integro-differential equations satisfied by the density matrices 
in this case can be obtained from (31) and (32) by making the summation 
for / to run from | to N. 


7. SPECIAL CASE 


In this section, we write down the explicit form of the correlated H.F. 
equations for the case of the simplest two-electron system, namely the Helium 
atom. We do not propose to evaluate the energy of the system as this has 
been done by several authors, but give the form of the equations (26) as this 
might enable one to have some insight into the meaning of these equations. 
We have here 


H.! = aE (x,) = @? {Pa - 2." 


2m ry, 


He? (%1, X2) = 07F (x, X2) + oB [F(r12), E (x) + E ()] 


4 20h Lew) + Bf (ra) [E (x:) + E (xa) 1S (rss) 


4 Bre e bs ¢ 12) (49) 


Te 
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Equations (28) become 
He! (x) #4 (%y) + [J He® (X15 X2) 2” 2) Axe] tr Or) 
— [J He® ra, 2) e* (%2) tr (Xz) Axe] Ye (0) 
= E,y, (1) (50) 


and a similar equation for electron 2 which may be obtained from the above 
by interchanging the indices | and 2. In the above we have chosen a dia- 
gonal representation for the matrix (Ajj), and E, and E, denote the diagonal 
elements. H,' denotes the kinetic energy operator of the clectron | plus 
its potential energy in the field of the nuclei. The second term in the brac- 
ket gives the average of the operator H,? (x;, x2) for all positions of the electron 
2 and the third term gives the exchange effects. Since H,? (x, x.) contains 
aterm a?F (x,, x2), the part of it containing the factor a* can be interpreted 
as the potential energy of electron 1 moving in the average electric field of 
electron 2. The terms af [/(r,.), E(x,) + E(x2)] give the influence of 
correlation on the kinetic energies of the electrons. The term [af f f(r,») 
bo” (X_) dx2] E(x.) especially shows that the kinetic energy of electron 1 
is not independent of the motion of the electron 2, but is correlated with it. 
Since correlation reduces the chances of two electrons coming close to each 
other, it also reduces the fluctuations in the kinetic energies of the electrons 
and thus tends to make the distributions of the kinetic energies rather uni- 
form. If we take f (ry) = ry, then af.* ff (1y2)/rie Ho? (Xo) dx, = aB-? and 
therefore the Hamiltonian for the electron | contains a constant term. By 
transferring this term to the right-hand side, the latter becomes («, — ae?) x, 
and thus the influence of the correlation factor is to reduce the energy of 
the system. Since f(r,.) is an increasing function of r,,, the effect of the 
term Be? f 2 (r,2)/ry2 in He? (x;, X2) will be to increase the average inter-electro- 
nic distance and consequently to diminish the energy of the system. 


8. NON-STATIONARY SYSTEMS 


There are several problems of physical and chemical importance that 
depend for their solution on the time-dependent Schroedinger equation 


(w+ A2)x=0 sy 


In this section, we consider such non-stationary systems and find out the 


orbitals (x,t) that make the function ® = g¥, where ¥ is given by 
A2 
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(3.a), the ‘ best’ approximation for the equation (51). These orbitals are 
given by the variational principle 


8 for (H+ 3) Od = 0 (52) 


8 f v* (He — iG 5.) Pdr = 0 (52a) 


G 
“+s ae Gi) +5 2; G (i, j, k) 


+ 4! Gi, j, k, l) 


ijkl 
G (i,j) = 2aBf (rij) + BF? (ray); 
G (i, j,k) = 287 [ F (rig) (fin) + rin) > fried) Sind} | 
G (i, j, k, D= 2B? {7 (rij) / (rea +S (rik) F (ru) +/ ‘(ria S (rie }. (54) 


We shall now evaluate the variation in (52 a) in two steps: first find out the 
variation of J, = {| ¥* H,.dz and secondly calculate the variation of 


Pb od 
7 f Gye dn. (55) 


(a) First consider 
5J, = f by* H,.¥dr + J 3YH.¥* dr. (56) 


Now f ¥* H. Wdz is given by (22) and f 5¥* H,.Wdr can be obtained 
from this by varying only the primed quantities (or the quantities with an 
aesterisk) without varying ¥,, %e, ...,%y. We shall denote variations of 
this type by 8*. Similarly {| 5¥?H,_'’* dz can be obtained by replacing the 
r’s in (22) by their complex conjugates and varying only the’ functions 
fy, Ya, ---5%my in I*(x,'xq ... x1’ | x... x) (= 1,2,....,6.. As men 
tioned earlier we follow the convention that the- operators He! and 2/2 
operate only on the unpeimed ae and that after these operations have 

N) and ¢’ =f in 
the integrals. 
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Since by (8 a) the I’’s are functions of p'j;, we have 


5* I’ (x,'x,’ eee xy’ | Xy eee x1) 


a i 5* 0! 
= Pat ce id 
oP ij 


4 


=i D> Pl; ij (> Soy * (x4’) xe «)) (57) 


{,j=1 





where p’1;ij is the co-factor of p’j; in p’1. 


Now 
[8¥*H.¥ dr 


= $* z tt... xP Gy... | %..- aa....am 
t=1 


8 i N ; 
= Ef Hel (a... x) ZF plrsag | Fb du* Ci) de (x;)| 


t=1 


xX dx, ... dx} 


z ng He Oe ss 
=P He 1 3 Boa" le 
+ z > P'1; aj 2 buby,* (xi') br ch (58) 


smi 


Now the terms corresponding to the cases i = 2, 3, ..., / can all be obtained 
from the first term of the above bracket by means of the cyclic permutations 
(x;, Xs); (X3, X3)3 --- (%, x1) of the integration variables and thus the sums 
corresponding to the cases i = 2,3, ...,/ are all identical with the sum cor- 
tesponding to the case i=1. Thus 


[ 8¥* HY dr 


= D nS He soo M0 
spp 2 P'1; 15 Sk (Xy') Ye a dx, ... ax 59) 


j=. 
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A similar expression follows for { 5¥H,%* dz in which the variation acts 


on the unasterisked quantities alone. Thus 


~— ScnS bap veieees 


x {He! (x1... 1) p's aye (xj) dx, .. . dx} 
4 [ SWH.P* dr. 


(b) Next we have 


* 
J, = f Gd ¥* oF dr — f cov ~ dr. 





Now 
5* f/f GY* = dr 
= a > f 8% * Gn, 1) =EEED ay, 
+ d | 6* nS GO vee xi) “EM de, — 
t=2 ; 
N—/—1 
+e pO Fee... rede .. 
where 
P11 Pi2 er ay eee ee Pils 
P21 Poe ee ee ee Palas 
ae: igen «SRR pn ene 
Pli MD a cieirelay oars Pllia 
OP Liat dP’ Lit, 
| U =n YT} 0 





N 
= a? ; f by * (2X, 6) oe ee t) dx, + te + ty 
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where f, and f, denote respectively the second and third terms in the above 
expression. We shall first consider 7. We have 


l= = a iT 5* f G (x; eee x1) °F tet dx, eee dX, 
- I G( ) 2 sep. d. 
<== i Xe KY ny; Ploy GX... AXI4y 
-3 ’ : G(x, ... 
i! dx, : —— 
2ly , Plan: ij > Sxby,* (xi') bh ah (65) 


i=1 
Now the integrals of the sums corresponding to the cases i = 2,3,...,/ 
are all identical with the integral for the case i= 1. Further let us inter- 
change the integration variables x, and x;,, in the sum corresponding to 
the case i=/+1. Then we get 


om Dag Use o> BD Gy... Sa] 


N 
>, k 
> 4 PY by,* (x,) a dx, esata X11 (66) 
where 
bk (X14 t) bk (Xe, t) tae bk (X141 t) 
ae ects Aa (67) 
P’ List, 1 P' Lis, 2 tae P' Lis, 41 
Next 


al a G(x... x) 
= z ot 3 Fa 
x > s Qh; ij ( pa Sox * (x4) be «)) 
+ + Detins (= Bue* (x14) PE yt (68) 


i=1 
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Now 


Qh; T+1, jf = Pl+s; 141, 3° (69) 


Interchange as before the variables x;,, and x, in the second term in the 
bracket of the above integral. Then we get 


— = Naf a Buhne (24) dy... deri 


x \ —— Dd) Qs; 13 Fk (Xj) 


+ G(x... Xua) Sim ii 4 tee | (70) 
(c) Thirdly the orbitals 4; (x) satisfy the orthonormality conditions 
f ba* (x) 44 (1) dx, = 0 (71) 
or 
Mei J Sib * (xy) Wi (1) dx, + Ari J oe * (1) 5x5 (1) dx, = 0 
(i, kK =.1, 2, ...,N) (72) 
where Ay; are arbitrary constants. 


Subtracting (72) from 5 (J, — iJ.) which is given by (60), (64), (66) and 
(70) and equating to zero the coefficients of the variations 5y,,* (x,), we get 


set 1 
{H! (x3) — dha? 2} dae (x) + 9 gay 
l=2 
xf Helo.) DS" pases dhe Op dite. de 


* am >? F f {IG (x, ... x) + G(x... xuv} 


pk 
aiid dX eee AX, 
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x f } (xX, ... x1) 2 Qh; 13 Ye (5) 


jal 


t 
r) 
+ G(x ..- Xt) » Pl4a; 43 “a ++ Ais 


j=1 


N 
= Z Ani Yi (%1) (73) 


and a similar equation for the complex conjugate of x (x). 


In the above equation, the Lagrangian multipliers A,; are functions of 
time; further when correlation is neglected, (21) reduces to the time-dependent 
H.F. equations as can be seen by putting a = 1 and B = 0. 
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SUMMARY 


The paper deals with a study of correlation effects in many-electron 
systems. Coulomb correlation is introduced into the theory by multiplying 
the Slater determinant formed from the one-electron orbitals by a correlation 
factor which is a symmetric and increasing function of the inter-electronic 
distances. The integro-differentia! equations satisfied by the best one- 
electron orbitals have then been deduced for non-stationary systems. From 
the extended Hart ece-Fock equations given by Léwdin, the integro-differ- 
ential equations satisfied by the density matrices have been derived. An 
expression for the energy-matrix of the system which is helpful in deriving 
a correlated Thomas-Fermi charge distribution, has also been given. 
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APPENDIX | 
We shall here evaluate 


7, - dp’ 
nt J Geax. x1) Sp OO. - Bey (1) 
and prove the relation (9) of Section 2. We have 


n= Q pot. x1) 


Gh. .+ 


j t N 
x l > S PN; ij +s +? ) | > | PN; ij =a (2) 
i=l j=1 tole = 
= T, + Tz (say) (3) 


where T, denotes the summation of i from 1 to 6 and T, denotes the sum 
of terms for i=/+1 to N. Since G(x, ..., x ) is a symmetric function 
of x;, ..., 1, the integrals of the sums corresponding to the terms i = 2, 3, 
...,/ are all identical with the sum for the case i= 1. Thus 


N N 
] P , . 
T, = (1) IG (x, ... x1) > PN; 19 ae dx, ... dXy. (4) 


Again the integrals corresponding to the cases j = /+ 2,/+ 3,...,N are 
all identical with the value of the integral] for the case j = / + 1, as can be 
seen by the transpositions of the integration variables (x14, Xis2)3 (Xti3» Xtis)s 
++, (%Ly, Xn). Further when j < (/ + 2), 


J pw; 3 dx1,2 ... dxxy = (N—D! Pl+1; 1j- (5) 
Thus 
T ~~ | IG (x, ... x1) 
iT dx, ... Ms 









0p’; dp’ 
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We shall next consider T,. 


By means of the transpositions (x1; X1;2) ... (X41, Xn) Of the integration 
variables, one can see that the integrals corresponding to the cases i= /+ 2 
to N are all identical with the value of the integral for i=/+ 1. Thus 


= ] , a Nfs. as P Lit, 5 , - 
Ts —_ Nt (N /) G(x, cee Xv) PN; lj yt dx, cee dxy 


(7) 
Dis 
l 
= ery ( N—) f GG... x0 
1+1 , 
0 ; , 
2 > Paztisg oe + OS — EY 
dp’ 2 
x P Nil 1, 14-2 =a | dx, eee dxx (8) 


which follows again by the same argument of permutation of integration 
variables used above. Integrating the first term in the above with respect 
to the variables x1,. to xy and the second term with respect to the variables 
Xu3 tO Xn one gets 


5 = l Gey .... 2) 
" /! dx, ... aXi4 


1+1 
d , i 
x jo ime l) , Plit; +159 ee + (N cae {— 1) 


/ P Lia li2 7. 
X Pl+2; l+1,1+-2 Sn dx, , . (9) 


Let us now consider the second term of the above expression and denote it 
by T,2. Expanding the determinant p},2;14:,1,2 along its last column, inte- 
grating with respect to the co-ordinate of the electron (/ + 1) and changing 
the integration variable x1. into x1, in the resulting integral, one gets 


T. = — Wa fF Gi. dda. deus 
A3 
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Q f l+ d tly b+ U 
x | Puss; iH eu a4 p _o : 7 p' ta1; Li, nat . (10) 


From (9), (10) and (6), we get on summing T, and T, and grouping together 
like terms that 


(T,+T) =) f IG (x, ... x1) 


1+1 
> y dp 
x e Plii; 15 — as) dx, eee AX 144 


_ os 
+ i dx, ... Xu 


«Sor l) 5 P Lia; Lit, j has 


’ dp 
+ Plat; Lia, Lia om (11) 


or 


| 
T; = i f GG, ... may, ... Gna 


«1D iit om Stas Pls 
a aP teed 12 
+ ¢ Plaa; l4nj (12) 


On writing the summations in the above term in a determinental form, we 
get the relation (9) of Section 2. 
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1. INTRODUCTION 


THE infra-red absorption spectra of the various isomers of the 0, m, p cresols, 
xylenes, ethyl toluenes and dichlorobenzenes have been investigated by 
Whiffen ef al/.! Their spectra reveal the existence of intense absorptions in 
the regions 700-800 cm.-' and 1200-1300cm.-' The region 700-800 cm! 
is characteristic of benzenoid substances, with the substituents attached to 
the nucleus. It has been found that there is a shift of this ‘aromatic band’ 
towards higher frequencies in passing along the series o-m-p in most of the 
disubstituted benzenes regardless of the nature of the substituents. 


The Raman spectra of the three isomeric cresols were investigated 
by Biswas,” in both liquid and solid states. The Raman spectra of o-p- 
chlorotoluenes and of o-p-xylenes in liquid and solid states were studied by 
Sanyal* and Biswas* respectively. According to Sanyal, p-chlorotoluene 
does not yield any line with a frequency below 306 cm.-! whereas there is 
an intense line at 163 cm.-! in the case of o-chlorotoluene. This line, he 
has attributed to the presence of dimers in o-chlorotoluene. 





Very little work has been done in the case of these isomers with reference 
( to the intensities of the Raman lines. A study of the frequencies and the 
intensities of the various Raman lines with reference to the position of the 
substitution has been attempted in the present investigation. 


2. EXPERIMENTAL 


F Extra pure samples of the substances distilled at constant boiling-points 
i have been used in the present investigation. A Hilger E 612 two-prism glass 
f 195 
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spectrograph, in conjunction with a Raman Source Unit, was used. The 
spectra of each series have been recorded under identical conditions. A set 
of graded intensity marks were also given on the same plate and _ the inten- 
sities were computed using a Moll microphotometer. 


3. RESULTS AND DISCUSSION 


Tables I to [V contain the frequencies and the integrated intensities of 
the Raman lines in different isomers. 


TABLE I 


Xylenes 





o-Xylene m-Xylene p-Xylene 





Frequency Intensity Frequency Intensity Frequency Intensity 
cm,—! cm.-! cm.-! 





170 
313 
459 
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TABLE II 
Cresols 
o-Cresol m-Cresol p-Cresol 
Frequency Intensity Frequency Intensity Frequency Intensity 
cm.-! cm. cm. 
190 3-266 
214 
274 1-431 0-476 337 3-130 
239 
445 2-946 455 0-468 464 1-980 
584 1-776 612 0-084 646 0-919 
735 3-012 
749 11-298 841 8-379 
833 1-931 
1002 2-190 
1044 4-407 1042 0-211 1100 0-227 
1081 0-629 
1252 3-390 +4 on 1251 2-083 
1378 1-209 1378 0-772 1380 1-740 
1400 2-625 1400 2-680 1400 2-638 
1598 1597 1598 
2-396 3-549 4-873 
1620 1614 1612 
2920 3-279 2929 2-292 2927 2-857 
3057 3-922 3057 5-746 3048 6-850 





From the results given in the above tables the following conclusions 
can be drawn: 


(1) The intensities of the corresponding lines are, in general, low in 
the meta-isomer. 


(2) In xylenes.—(a) The strong line in the region 700-800 cm.—! repre- 
sents the benzenoid vibration. According to Whiffen and Thompson this 
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o-Chloroaniline 





Frequency 


cm.-? 





168 
257 
467 
555 
678 
833 
1020 
1304 





TABLE III 
Chloroanilines 
Intensity Frequency 

cm.~! 
1-900 193 
1-473 244 
1 -433 408 
1-008 530 
1-415 694 
1-878 893 
1-843 1073 
2-231 1300 

TABLE [V 
Chlorotoluenes 


o-Chlorotoluene 





Frequency 
cm. 





164 
247 
675 
801 
1040 
1204 
1590 
2920 
3062 


Intensity 


1-692 
2-383 
1-167 
1 +366 
1 +243 
0-947 
1-545 
1-149 
1-543 





m-Chloroaniline 





Intensity 





tu 


-060 
| +362 
| +366 
0-683 
1-641 
1-071 
| -092 
4-271 





p-Chlorotoluene 











Frequency Intensity 

cm.~? 

307 1-180 

634 0-801 

796 1-803 
1090 | -866 
1208 0-639 
1608 9-819 
2927 0-497 
3050 1-078 
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line in the infra-red appears at 743, 765 and 790 cm. in passing from ortho- 
to the para-compound. In the present investigation strong lines have been 
observed at 733 cm.-!, 724cm.-' and 784cm.— in the three isomers respec- 
tively, thus showing a diminution in frequency in the meta-isomer. 


(b) The line in the region 2880-2860 cm.-! represents the symmetrical 
stretching of the methyl group attached to the benzene nucleus. This line 
has a2 maximum intensity in para, but it is very weak in meta and is absent 
in ortho-isomer. 

(c) The line at 2916cm.-!, representing the asymmetrical stretching 
of the methyl group, shows a gradual increase in intensity from ortho to 
para-isomer. 

(3) In cresols—(a) Strong lines at 749, 735 and 841 cm.-! have been 
observed corresponding to the benzenoid vibration, while Whiffen and 
Thompson recorded it at 752, 775 and 827 cm. in the o-m-p cresols res- 
pectively. 

(b) The intensities of the C-C vibration at 1000cm. and the CH; 
asymmetric stretching in methyl attached benzene ring at 2916cm.— behave 
in the same way. 


(c) No line has been observed in the region 2880-2860cm.-! which 
represents the CH, symmetric stretching mode. 


(4) In chloroanilines.—{a) There is a slight increase in intensity of the 
C-Cl vibration at 678cm.-! in the meta-compound. This observation is 
different from what has been found in the other substances. 


(b) The line at 1300cm.— in meta-isomer shows a doublet structure 
as a result of which the integrated intensity has been considerably enhanced. 


(5) In chlorotoluenes the intensities of the lines in ortho-compound 
are greater than in the corresponding para-isomer. 


The para-chloroaniline and meta-chlorotoluene were not investigated 
due to non-availability of the substances. 


Since the vibrational intensity of a Raman band is a dependent func- 
tion of the derivative of the polarisability, which is a related function of the 
electronic state of the bond, the observed variation in intensity in the three 
isomers may be explained, by considering the inductive effects brought into 
play by the various substituents attached to the benzene ring. 


Shorygin® while studying the influence of the second substituent on the 
mono-derivatives of benzene has observed that the intensity of the first sub- 
Stituent depends upon the electronegativity of the second substituent group. 
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In the case of the intensity of the nitrogroup band, the intensity increases 
very sharply when electropositive radicals are introduced in the para posi- 
tion. He has also pointed out that the halogen atoms behave as electro- 
positive groups as far as the intensity of the Raman scattering is concerned, 


The compounds studied in the present investigation have the following 
groups: CH, groups in xylene; OH and CH, in cresols; NH, and Cl in 
chloroaniline; and CH, and Cl in chlorotoluene. 
known to be electropositive. 
for each group. 


All these groups are 
The degree of electronegativity is different 


When an electropositive group like CHs is attached to the benzene ring, 
this group will force the electron into the ring due to the electron repelling 
character and consequently increases the electron density at the ortho and 
the para positions with respect to it and the effect at the meta is much less. 
If another methyl group is substituted in the para position, as in para-xylene, 
the positions a and b (Fig. 1) possess high electron density owing to the 
repulsive inductive effect of the alkyl group in the ortho position. Thus 
both the groups aid each other and the electron density will be increased. 


CH, 
! 


LN 





CHs 
Fic. 1. 


When the methyl group is attached to the meta position, due to cumula- 


tive effect of the electric charges induced in the ortho and the para positions, 
the electron density will be maximum. 


In the ortho position, the inductive effects due to both the groups oppose 


each other and consequently change in the electron density induced in the 
ring will be minimum. 


Thus it is seen that the changes in electron density induced in the ring 


follow the order meta > para > ortho. This will be the case with all the 


electropositive substituents though the changes may vary depending upon 


the degree to which they are electropositive. The observed results may 


thus qualitatively be explained on the basis that the increase in the electron 
density may cause a decrease in intensity of the Raman lines. 
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4. SUMMARY 


The Raman spectra of the o-m-p isomers of some substances have been 
investigated. It has been found that in all the cases the intensities of the 
corresponding Raman lines have a maximum value in the ortho-compound 
and a minimum value in the meta-isomer. An attempt has been made to 
explain the observed results on the basis of the change in electron density 
caused by the inductive effect of the various substituents attached to the 
benzene nucleus. 


One of the authors (K. R.) is thankful to the Council of Scientific and 
Industrial Research, Government of India, for the award of a Junior Research 
Fellowship. 
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1. INTRODUCTION 


THE position of C = N stretching vibration in nitriles has been extensively 
studied by various investigators,'~* and it is found that the above character- 
istic vibration occurs in the region of 2230cm.~! 


The infra-red absorption study of solvent effect in relation to inter- 
molecular forces has been made recently by Caldow etal.4> They have 
pointed out that the variations in the vibrational spectrum may reveal altera- 
tions in the bond character caused by the particular type of intermolecular 
interaction and thus may even suggest a method of union between interacting 
partners. ; 


A given vibrational chromophore may also show different behaviour 


according to the nature of the rest of the molecule in which it occurs. Jesson 
and Thompson® observed in the infra-red that on passing from vapour state 
to solutions in CCl, and CHCl,, the C =N group vibrational band of 
CH;CN steadily increases in intensity whereas that of the C = N group in 
CCI,N steadily decreases. It has also been noticed that in these more com- 
plex cyanides C = N group frequency is almost independent of the solvent 
although the band intensity varies considerably. 


It has been observed by Caldow ef al. and Thompson and Steel’ that 
the infra-red intensity of the C = N vibration of benzonitrile in CHCl, is 
approximately double that of the same in CCl,. A study of the solvent 
effect on the infra-red intensity of the C =N stretching vibration for some 
nitriles made by Foffani et a/.* shows that the band intensities for the nitriles 
increase with the polarity of the solvent except in ketonitrile where the trend 
is reversed. They have attributed this difference in behaviour to the pre- 
sence of two competing bond moments in ketonitriles. 


It is clear from the above that, although some work has been done in 
infra-red, no work has been done regarding the study of the solvent effect 
on C=N vibration in Raman effect. Accordingly, a study of the effect 
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of various types of solvents on this characteristic vibration of nitriles has 
been taken up in this investigation. 


2. EXPERIMENTAL 


All the liquids studied in the present investigation have been distilled 
at constant boiling-point before use. The Raman spectra of various mix- 
tures (50°{ by volume concentration) have been photographed under identical 
conditions using Hilger Raman Source Unit along with E 612 Spectrograph. 
A set of intensity marks are also recorded on the same plate and the inten- 
sities estimated using the Moll recording type microphotometer. 


3. RESULTS AND DISCUSSION 


The ratio of the integrated intensity of the C = N frequency of benzo 
and acetonitriles in various solvents and the 458 cm.~' line of pure CCl,, 
after applying the refractive index correction suggested by Woodward and 
George,® is calculated and the values thus obtained are presented in 
Table I. 


TABLE I 
Solute Benzonitrile Acetonitrile 
a 2224 cm.-! 2250 cm.-} 
Solvent (50% by 
vol. conc.) 
Pure es 2-320 0-238 
Carbon tetrachloride si 1-172 0-092 
Chloroform 0-940 0-078 
Bromoform 0-172 0-103 
Ethyl acetate 0-916 0-099 
Methy! butyrate 0-512 0-106 
Ethy! butyrate 0-913 0-100 
Benzene 1-472 0-120 
Toluene = 1-349 0-122 
Ethyl benzene on 1-482 0-101 
Chlorobenzene re 1-425 0-133 
Bromobenzene me 1-318 0-104 
Dioxan ‘a 1-269 0-096 
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From the results given in Table I, the following conclusions may be 
drawn: 


(1) In the case of benzonitrile mixtures, the intensity of C =N vibra- 
tion does not show any appreciable variation in the aliphatic non-polar 
solvent CCl, whereas it is considerably d'minished in the aliphatic polar 
solvents. It shows an increase in intensity in the case of all the aromatic 
solvents studied. 


(2) In the case of acetonitrile mixtures, the intensity of C = N vibra- 
tion does not show any variation in aromatic non-polar solvents like benzene 
and toluene. In both polar and non-polar aliphatic solvents and in polar 
aromatic solvents, in general, an appreciable diminution in intensity has 
been obtained. 


(3) In heterocyclic non-polar solvent dioxan an increase in the intensity 
is noticed in benzonitrile mixtures while there is a decrease in intensity in 
acetonitrile mixtures. 


Thus it may be inferred that the aliphatic and aromatic non-polar sol- 
vents have no effect on the aromatic and aliphatic nitriles respectively. In 
all other cases, the solvents may interact with the solutes thereby altering 
the bond length of the chromophoric group and hence affecting the polar- 


isability derivative. This may probably be the reason for the variation of 
intensity noticed in those cases. 


According to Caldow etal., the C=O and C=N groups should 
behave roughly in a similar fashion as regards interactions with solvent 
molecules. The observations made by Venkateswarlu and Ramaswamy" 
show that the ratio of the intensities of C = O vibration in CCl, and CHCl, 
solvents in general is found to be approximately 1: 2 which is in good agree- 
ment with those made by Caldow eta/. But in the present investigation 
the behaviour of C = N vibration is found to be different from that of the 
same in the infra-red in the sense that the intensity of the line in CCl, is 
observed to be great when compared with that in CHCl, the ratio of the 
intensities being roughly 5: 4. 


4. SUMMARY 


The effect of various solvents on the intensity of C = N vibration of 
benzo and acetonitriles has been studied. It bas been observed that the 
aliphatic and aromatic non-polar solvents have no effect on the aromatic 
and aliphatic nitriles respectively. In all other cases, a marked variation 
in intensity has been noticed. The intensity of C = N vibration in Raman 
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effect shows a diminution when the solvent is changed from carbon tetra- 
chloride to chloroform while it shows a considerable increase in infra-red. 


One of the authors (C. B.) is thankful to the Government of India for 
the award of Senior Research Scholarship. 
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ABSTRACT 


It is shown that the interaction between two charged particles 
reduces to that of a potential with one particle when the mass of the 
other tends to infinity. Recoil corrections when one of the masses is 
finite and large is also estimated in the case of two fermion scattering and 
bremsstrahlung in two-fermion interaction. 


1. INTRODUCTION 


In classical mechanics it is well known that when two particles collide there 
is only a change of momentum without transfer of energy if one of the 
colliding particles has infinite mass. In electrodynamics the interaction 
of a charged particle with a potential is characterised by a change in momen- 
tum without any change in energy if it is assumed that the seat of the poten- 
tial is stationary during the process. However if the potential were due 
to a particle of finite mass its motion will have to be taken into considera- 
tion and this implies an exchange of energy also. For example, this is the 
case in electron-electron scattering where the interaction occurs through 
the exchange of a virtual photon or equivalently through an instantaneous 
Coulomb potential and transverse photons. Hence it is found interesting 
to investigate whether in the case of a two-particle interaction when one of 
the masses tends to infinity the situation will reduce to that of an interaction 
with a potential. This is indeed found to be so in correspondence with the 
classical case of change of momentum without transfer of energy. If 
instead of assuming the mass of one of the particles to be infinite, we con- 
sider it to be large compared to the other, then we can obtain a recoil 
correction proportional to the ratio of the masses. 


* Permanent Address: University of Madras, Madras, India. 


t Council of Scientific and Industrial Research, Government of India, Junior Research 
Fellow. 


§ University Grants Commission (India), Junior Research Fellow. 


306 





A Limiting Process in Quantum Electrodynamics 207 


2. SINGLE PHOTON EXCHANGE BETWEEN TWO FERMIONS 
We shall first demonstrate that a two-fermion interaction reduces to 
a potential scattering of a fermion when the mass of the other tends to 
infinity. The Feynman diagram for the interaction of two fermions of 
masses m, and mz is given in Fig. | 


ae 


Fic. | 


The matrix element for this process is 


M = 4ne? “ pe [@ (Ps) yuu (p1)] a [4 (ps) vult(pa)] 


Spins 1, 2, 3, 4 
with 
—Pi =k = )2—Fa (1) 
which is the conservation law obtained from the space-time integration. 
From purely kinematical considerations, making use of these conservation 
laws it can be shown that if m, tends to infinity and the initial momentum 


P, is zero, then the fourth component of the final four momentum p, is the 
same as that of p,, i.e., the rest mass. Hence the fourth component of the 
four momentum of the virtual photon is zero. However, taking m, to be 


finite and large compared to m, we find that D, does not vanish and hence 
there is a recoil. This recoil correction is expressed conveniently in terms 
of the ratio of the masses which vanishes when m, tends to infinity. The 
square of the matrix element (1) reduces to 


4re? Pitm_ &s+m 
| 2 1 1 3 
amd -|% ae [as 2 5 8? [re ge se] 


1 Pot+m. Pat my, 
oe Nd [% 7B, 7s ral 


here is no recoil. 
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Hence the cross-section becomes 


o = transition probability per second x I/v 
= | 2m | M |% (E) dE (3) 


where v is the velocity of the particle (= p/E) and p(E) dE the density of 
final states of the single particle. 


Since E,= E,=E and |p.| =|p,|=|p| and if @ is the angle 


between p. and p, we have 


4etE? sin?6 
o= Sy [1-8 >|: (4) 


This is identical with the scattering cross-section of a particle by a potential. 
Thus we find that the second order process in a two-particle case reduces 
to a single scattering by a potential when the mass of one of the charged 
particles tends to infinity. The cross-section with recoil is given by 


2 
Gr ecoil an —%) + - = * Ee x {! * E, (1 _ IES 


(i) 


< (Et) + a + ate) 9 
and 


— {sm [E.+ Z + E, cos 0] 


cos@  2E, cos @ m,° 
—m(g, + Sp, + “ae ) — ogee} 
where « = m,/m, and terms of the order 1/m,? have been neglected. 


3. Two-PHOTON EXCHANGE BETWEEN Two FERMIONS? 


We now consider the process of higher order where two photons are 
exchanged between two fermions. The two possible Feynman diagrams 
are given in Fig. 2. 


The matrix element can be written down as usual. But in this case an 
integration over the four momentum k of the photon also occurs. This 
is due to the fact that there are two virtual photons exchanged and the 
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given transfer of momentum between the two fermions can be shared in 
any ratio between the two photons. The net transfer of energy between the 
two fermions is zero in the limit m, tending to infinity and this gives rise 
to the factor 5 (q4) in the matrix element. Thus the matrix element for the 
R diagram is 


Mem[ OD 


B, Spins 1, 2,3,4 


a Pitk tm atk 
x { f [#0 fy zP— ma Yat (Ps) | 2 Og — Ke 


x E ()4) 'y is aa Be. 42 ye Ya (Ps |}. (6) 





This integral is divergent but can be evaluated by the usual tec‘inique* of 
introducing a small negative imaginary part ie and making ¢ tend to zero 


and also the convergence factor A?,,,.. Thus 
Eps A. 
e+ * S dk 
“7 a) a (Pa) =e — m,2"'#" (P| 2k 
i dk a 
rs (k coe Was \(q 7 k)? a Tail ; “ 


It is to be noted that g is only a three vector in this expression (7) since the 
conservation laws give rise to the factor 4 (q,). 
Ad 
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The corresponding Feynman diagrams for the double scattering of a 
fermion due to a potential are given in Fig. 3. 


oa ww) 
/ bee \ 


ws — i 


ila mas 
V(k) ig ; Ma 
[s} A 

[x] ~. ; he 


Fic. 3a Fig. 36 


Ay 





The matrix element for the R’ diagram is 


OOD f [evan Seek tm e-wtoa] / 


Spins 2, 4 


x ae 
(%2— Ma.) (@ — 0? — Mane) 


Here the integration is only three dimensional because the conservation 
laws yield the factor 5(k,)8(qg4— k4). It can be shown that expressions 
(8) and (7) are identical by carrying out the integration over k, in (7). This 
is performed by taking the principal value of the integral which is equal in 
magnitude and opposite in sign to that obtained by performing the inte- 
gration along a semicircle of radius « around to pole ky =0 as «—0, 
provided there are no poles in the upper half plane as is the case here. Thus 
we have 





(8) 


| ¢ dk, 
“a2 “"? 





Hence the expressions (8) and (7) are identical. Similarly it can be shown 
that Ms and Me become identical in the limit m,—> co and thus two-photon 
exchange is found to be equivalent to double scattering by a potential in 
this limit. 
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4. BREMSSTRAHLUNG IN TWO-FERMION INTERACTION 


In this process we shall demonstrate that bremsstrahlung becomes more 
important as the mass of one the particle increases and in the limit m,—> ©, 
the process reduces to a pure bremsstrahlung. The Feynman diagrams 
are shown in Fig. 4. 





Fyc. 4a Fic. 4 b 


The matrix element for the R diagram is given by: 


ame A (Ti (ps) yu (Prd) 


Spins 1, 2, 


oe [#(p0) zg (p -. 3 jt m,2 Yult (Pa) } . (9) 


” The square of the matrix element is then given by 


ity = [4 ETD»  aeee/ 


Cree ee 


xy Pat ms. .+K +m, g (@ 2) 
oe  e "te, + Kk) — m* 2E, 


J Sp [-. Pi, +m. Bs ]} 


“9 7E, ‘* 25, 
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[Sr] ; Pat K +My, P2tm, 
tier + #)?— m2" ~2E.- 

5 Pot & +m, Pat me, 5 
XTep se ome” 3B, I['+ aml} 7 


where 6 is the increase in energy of the mass m,. In the limit m,—> oo, the 
factor 5(k,) can be obtained as in the case of single photon exchange and 


> 
the four vector k reduces to k. 


The Feynman diagrams for ordinary bremsstrahlung process are given 


in Fig. 5. 
A 
va) ae 





(| 
fs] 
Fic. 5a Fic. 5b 


The matrix element of R’ diagram is 


_ 4re? i Fat KF +m, 
Ms -["]] dae eft A ‘wa. (11) 


Or, the square of the matrix element is 


2 
4re? q Pi,+k+m 
Ma 2 
paleo 


x y, fat mS Pot hk +m , 2 Fat). 








2E, “(Pi,+ kP— my? 2E, 
(12) 
On carrying out the spur calculations it is found that the expressions (10) 
and (12). Similarly the limiting forms of | M, |? and | My |? will turn 
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out to be identical. Thus we find that bremsstrahlung in a two-fermion in 
teraction reduces exactly to ordinary bremsstrahlung in this limiting case. 


The cross-section taking into account the recoil, i.e., when m, is large 
and finite, is given by 


> 2E 4 


[e+ 2p,.%]  (%—2,.%) 
2 


1 | e dds Pa-: , XE 
+5 |= | e? 4 ‘sin 0,d0, sin 0,d0,dp ——* 
an >, De 2d05 ,48,d¢d 5) 








Frecon 


— 2E,"E, — 2m,"E, + 2m,"E, 
> - >\2 
(ke +2p,.k ) 
2E,"E, + 2m,*E, — 2m,*E, 
+ > - >\2 
(x —2p,.k ) 








> [m,"E, — m,*E, — E,E, (E, + E,)] 

(2+ 2p,.k) (k®— 2p,. k) 
where ops, is the normal bremsstrahlung cross-section and X, the expression 
defined previously. 


+ (13) 


This work was done during the stay of one of us (A. R.) at the University 
of Bern under the auspices of Swiss Atomic Energy Commission at the kind 
invitation of Professor Andre Mercier. He wishes to express his gratitude 
to Professor Andre Mercier for his hospitality and the facilities afforded 
for work. G.B. is indebted to the Council of Scientific and Industrial 
Research and S. I. to the University Grants Commission for the award of 
Junior Research Fellowship. 


1 Throughout this paper we use the following notation. We denote the four vector p by p 
> 
and the three vector by p. The daggared operator # is defined as P = 1:Pt~ VePe— YyPyu~ VaPe 
> > 
and the scalar product of two four vectors as p, . fg = Py; Po: — Py » Pg» We have also used the 
natural units # =c = 1. 


2 We are grateful to Mr. A. P. Balachandran for suggesting this problem and drawing our 
attention to the fact that in spite of the crossing of photons in the four vertices the interaction 
in the limiting case should reduce to a second order scattering of a potenti]. 


3 Feynman,R.P., Phys. Rev., 1949, 76, 769. 
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ABSTRACT 


An expression for the charge density of proton is presented and the 
form factor calculated on the basis of this is compared with the experi- 
mental results. The results are also compared with those of the exponen- 
tial model and the cut-off model of Hofstadter (1959). 


INTRODUCTION 


THE experiments on the scattering of electrons by protons have yielded 
valuable information on the structure of the proton.!. The experimental 
results are interpreted in terms of the Rosenbluth expression? for the scatter- 
ing cross-section, modified by the introduction of the charge and the 
magnetic moment form factors F, and F,. 





do =f i gat 3 ° 1 
dQ ~ \JE Pre 2! 
0” sin! 5 + Face ° sin? 


x {F24 2S [2 2 (F, + KF,)? tan? 


, 0 
( 4M?2c2 2 


+ K*F*]} (0) 
where Ep is the incident electron energy, M is the mass of the proton, g = 
(2/4) sin /2 is the 4-momentum transfer and K is the anomalous magnetic 
moment of the proton. This formula is expected to be valid in the range 
of the electron energies and momentum transfers investigated at present. 
The phenomenological analysis of the data available is consistent with the 
simple assumption F, = F, = F, and that F may be interpreted in terms of 
a static charge distribution p(r) as 


mn * f p(r) sin gr . rdr. (2) 
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Several different charge distribution models may be written down, all 
of which give equally statisfactory fit to the experimental data. One of 
the best known is the exponential model which gives for the root-mean- 
square radius of the proton 0-80 x 10-*cm. All other satisfactory models 
also give similar values. Hofstadter* has pointed out that although the 
simplest meson-theoretic model 


e7ar 
p(r) = (3) 
is inadequate to fit the experimental data, this distribution modified for 
small values of r by the introduction of an additional cut-off parameter can 
satisfactorily fit the data. Hofstadter’s model is given by 


p = prt O<r<d 

—br (4) 
p= on d<r<oo 

and contains an unpleasant discontinuity at r=d. In view of the theore- 
tical interest of such models we have explored further in an attempt to write 
down a satisfactory well-behaved charge distribution function. Such a 
function is described in §2, and the results are discussed in §3. 


2. RESULTS 


The general aim is to change the expression (3) to 


pr) =[" ] 10 (5) 


r2 


where f(r) has suitable analytic properties. We require f(r)—»>0 for r—0, 
and f(r)—> | as r—» oo, so that the divergence of p(r) at r= 0 is removed, 
and piling up of charge at the centre of the proton is avoided. Our first 
attempt was to choose 


f(r) = (1 — e®°). (6) 


[t was however found that no choice of parameters a and f could produce 
a satisfactory fit to the experimental data. 


Next, we have tried a modifying function of the type 
f(r) = (1 — e 8")? (7) 


as we wish to keep the charge-density positive at all values of r. This com- 
bined with (5) and (2) gives 
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n 
n+ | 


F (q) = [zc + n)(2 +n) j tan? 4 — 2 tant ( 


i eal (. + 2 if 


where n = a/B. 


Now F? versus q? is plotted for different values of a and 8 varied systemati- 
cally. It is found that a satisfactory fit to the experimental data is obtained 
for 


a = B = 2-557 x 10% cm 
a“! = 0-39 x 10-3 cm. 


With these values of the parameters, one easily finds for the root-mean- 
square radius 


a= 0-85 x 10-'¥ cm. 


Figure | shows the comparison of the experimental values of F? and the 
values calculated on the above model. The errors on the experimental 








Fic. 1. Plot of F* vsg*. The abscissa is in the units of 10°° cm-*. 


values are ~ 10%. It can be seen that the agreement is quite good up to 
the highest momentum transfers yet observed, and a breakdown of the stati¢ 
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model is not yet in sight. The shape of the charge distribution p (r) multi- 
plied by r? is plotted against rin Fig. 2. Also shown r in this figure is the 


exponential charge distribution. 
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a 
The abscissa is in the units of 10-¥*cm. The solid line 


Fic. 2. Plot of 4wr? p (r) vs. r. 
repiesents the charge distribution medel given by equations (5) and (7), and the dotted line gives 


the exponential model. 
DISCUSSION 


We shall now make a few remarks on the results described above. The 
values of a and a confirm the general trend of the results of Hofstadter. 
The root-mean-square radius is quite close to that obtained by other models, 
viz., 0-80 x 10-1%cm. for the exponential model, and 0-82 x 10-14% cm. 
for the models given by equation (4). The value of the range parameter a~ 
is larger than that obtained in the exponential model, but is comparable 
to the value b-1 = 0-4 x 10-1%cm. obtained by Hofstadter for the models 
of equation (4). It is also clear that whatever the type of the modification 
introduced in the simple meson-theoretic charge distribution (3) for small 
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r it has to be operative for r<0-4 x 10-'* cm., as indicated by the vak 
d ~ 0:4 x 10-" cm. found by Hofstadter and a = B obtained by us. 
such distances, one would expect the simple z-meson cloud effects to bg 
modified by the effects of the K-mesons (4/M;,c ~0°4 x 10-% cm.) and 
the nucleon core (4/Mc = 0-21 x 10-'*cm.) structure. However, in view 
of the fact that the simple meson-theoretic model (3) would only be vali 
for rz 1x 10-cm. and the purely phenomenological approach of oug 
model, it would be unwise to stress the discrepancy between the ranges 
parameters obtained here (0-4 x 10-'%cm.) and required by the meson 
theory (0-7 x 10-*%cm.). § 


In view of the result a = 8, the expressions obtained from (6) and (7}j 
would appear to be power series in exp.(— ar). However, in view of our 
observation that (6) is quite unsatisfactory, one is tempted to interpret (TF 
as indicating a modification of the simple meson-theoretic potential! 
 ~exp.(—ar)/r to & ~exp. (— ar/2) [1 — exp.(—ar)]/r resulting frome 
the effects of the proton structure, 3 


Finally, we note that the charge distribution model obtained from i 
equations (5) and (7) [but not that obtained from equations (5) and (6)] rede 
ces to the exponential form exp. (— ar) for small values of r. : 
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